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I. INTRODUCTION 


Every quantum system unavoidably interacts with the surrounding environment. 
The dynamics describing such open quantum systems are in general very difficult 
to treat as they involve the whole history of the process. In order to analyze the 
physical features of these systems one needs to consider effective descriptions, which 
have a simpler structure and are then more suitable for a detailed analysis. If the 
timescale of the system is much smaller than the one of the environment, one can 
use the Markovian approximation. Markovian dynamics are the most widely used, 
and their mathematical and physical features very well knowiA However, many 
physical systems in chemistry and biology, do not display such a clear separation 
of timescales and they cannot be described by Markovian dynamics^—. In order 
to study such non-Markovian systems, one then needs to look for different kind of 
effective descriptions, not based one a timescale selection. The main problem one 
has to face when dealing with open systems is that the environment has a very large 
number of degrees of freedom which are hard to control. A possible solution to 
this problem is to consider an effective description which involves a small number 
of degrees of freedom. In this paper we prove that one can effectively describe the 
influence of a big environment, with an effective environment made of a (much) 
smaller number of constituents. 

The independent oscillators (10) model, is one of the most widely used to de¬ 
scribe open quantum systems^ 7 . In this model, the environment is described by a 
set of independent harmonic oscillators which are linearly coupled to the relevant 
system. The 10 model has been widely studied, and it allowed for the description 
of Markovian and non-Markovian quantum Brownian motion^ - -. Furthermore, one 
can derive a Generalized Langevin Equation, which is suitable for a phenomenolog¬ 
ical description of thermal and diffusive properties^ - —. However, though the GLE 
shows the existence of non-Markovian effects, it does not allow for an effective de¬ 
scription of these. In a recent paper—, it has been shown that a chain model for 
the environment suits better for this goal. In this model, which is widely used also 
in other helds^ 1 ^, the environment is described by a chain of Erst neighbor inter¬ 
acting oscillators. The ordering in the interaction (absent in the 10 model), gives a 
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clearer physical picture of how the non-Markovian effects build up and allows for an 
effective description of the dynamics. 

A first attempt of effectively describing the short time dynamics of a system has 
been made in^. Defining some “collective modes”, the authors obtain an effective 
Hamiltonian and they show that the cumulant expansion of this is equivalent to 
that of the original Hamiltonian up to the third term. Though they can argue the 
equivalence of the two dynamics for “short times”, this result is not rigorously proven, 
and it does not give the error made by considering the effective dynamics. 

A step forward the derivation of an effective dynamics has been taken inr^, where 
the authors solved the Heisenberg equations of motion for the chain models obtaining 
a description equivalent to the GLE. Unlike the GLE, which cannot be analytically 
solved for any environment, the dynamics obtained in-^ is completely general and 
its structure paves the way for an effective description. The aim of this paper is to 
complete this program: we prove that, for a harmonic system, one can derive an 
effective dynamics for any set of chain parameters, giving an explicit bound on the 
error one makes considering the effective dynamics instead of the full one. 

The paper is organized as follows: in Section II we summarize the main result 
of—, giving a mathematical account of it; in Section III we prove that one can give 
an effective description of the full dynamics, we prove that the error one makes is 
bounded. We eventually provide the relationship among timescale, number of modes 
and error. 

II. NON-MARKOVIAN DYNAMICS OF AN HARMONIC 
OSCILLATOR 

Our starting model is a system interacting linearly with an environment of N 
independent oscillators (10). Such a model is described by the following Hamiltonian: 



( 1 ) 


fc=i 


k =1 


where at this stage V(x ) is a generic potential, while later we will focus on harmonic 
systems, i.e. purely quadratic potentials, which one may see, if one wishes to do so, 
as Taylor expansion of the generic potential. Since we are considering only short 
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times such a Taylor expansion might be justifiable, but we do not do that in this 
paper. x,p are the position and momentum operators of the relevant system, and 
qk,Pk are position and momentum operators of the environmental oscillators. Such 
oscillators have proper frequency ojk and they are coupled to the system via the 
positive constants ty. 

The Heisenberg equations of motion for the environmental oscillators are easily 
written as follows: 

-^q(t) = -u> ■ q {t) (2) 

where q 3 = (qi, ■ ■ -q n ) is the vector of the environmental position operators, and 
c j is the diagonal matrix of the oscillators frequencies: u: = diag(a;^,... cwith 


(jJ\ < U>2 < ■ ■ ■ U>N- 


As we have already mentioned, a chain model for the bath is more suitable to 
study the short time behavior of the dynamics. Such a model is described by the 
following Hamiltonian: 

2 N N 

+ CM + DxX t + V D h I h X„ + £ - ( P * + a l X l) (3) 

k =2 k =1 


The chain oscillators have position Xk, momentum P*,, proper frequency f4, an d 
display a first-neighbor interaction with positive coupling Dk■ The Heisenberg equa¬ 
tions for Xk read 


dt 2 


■X(t) 


—T ■ X(t) 


( 4 ) 


where X 3 = (Xi,... Xjv), and T is the following tridiagonal matrix: 

/ O? -Di 0 ... \ 


T = 


—D 1 -d 2 ... 

o -d 2 


( 5 ) 


V : : 

In order for the chain description to be equivalent to the initial 10 model, one needs 
to require that the dynamics given by (pQ) and (J3D are equivalent. Accordingly, the 
parameters entering (J3J) are not free, but particular combinations of the parame¬ 
ters of and the chain oscillators A" must be specific linear combination of the 
independent q. We define 


Xj — Ojkqk ■ 
k 


( 6 ) 
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where O is an orthogonal N x N matrix. The condition of equivalent dynamics is 
fulfilled only when the systems of equations of motions fl21) and (JH) are equivalent. 
Substituting (l6|) . one finds that this is true when 


T = O uO T 


(7) 


The problem of determining a matrix starting from its eigenvalues is known with the 
name of Inverse Eigenvalue Problem!*^. For our case of study, such problem can 
be exactly solved, i.e. one can analytically determine the entries of O, and through 


them the matrix T and the parameters of the chain 0*., D For the details of the 


solution of the IEP and the derivation of the parameters one can refer to^, here we 
summarize the result. The entries of the matrix O read 




where Pj( A) is the characteristic polynomial of the j-th leading principal minor of 
T, evaluated in A. The explicit expressions for the Pj are determined recursively 
exploiting the following recurrence relation: 


fWA) = (&j - A)Pj(A) - DjP,_,(A) 


(9) 


with P—\ = 0. Once the transformation matrix O is determined, the entries of T 
are given by the following relations: 



( 10 ) 


k 



( 11 ) 


k 


From now on we will consider the matrix T as known, i.e. as fully determined in 
terms of the parameters of the 10 model, and we will assume the dynamics given by 
Eq. (J3]) to be equivalent to that given by Eq. ([I]). 

A. Non-Markovian dynamics of an harmonic oscillator 

After having established the equivalence between the two models, we determine 
how the chain modes affect the dynamics of the relevant system. At this purpose, 
we need to solve the set of equations of motion (J4]). Since the equations of motion 
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are equivalent for a quantum or classical systems, we treat the latter case in order 
to keep the treatment easier. We consider the relevant system to be an harmonic 
oscillator with proper frequency hi, in such a way that (J4j) can be explicitly rewritten 
as follows: 

A 2 

—x(f) = -tt 2 x(t) + DXi(t) (12) 

^W(t) = -nfXdt) + + DiXi+i(t) , 1 < i < N - 1 (13) 

-j^X N (t) = —Q 2 n X N ( t) + D N _ 1 X N -i(t) , (14) 

We independently solve the equations of the previous system, and rewrite them as 
follows: 

x(t) = /o + f Sm[n f~ s)] DX 1 (s)ds (15) 

Jo “0 

x,(t) = fM+t — 7 ) ~ 1 < i < N—l ( 16 ) 

Jo 

Av(^) = /n + f —^-—Av-iAv-i(s)ds , (17) 

Jo “JV 

where we have relabeled hlo = A Xq = x, and 

/i(t) = Xi(0) cospirf] + W(0)^&fl . (is) 

In order to obtain an equation for x in terms of the X i: we substitute recursively 
in x(t) the equations for X\ (t), X 2 (t),... X n (t), and we prove the following Theorem. 


Theorem 1: Let {x{t), Xi(t)}i = i^__ : N 6 R be the set of functions solving the sys¬ 
tem fl 1 5H - fl 1 7[) . Let the functions K^t — s) : M 2 —> K. and f t (t) : M 2 —> R. be dehned 
recursively as follows VI < i < N: 

Ki(t — s) = [ Ki_i(t — l) sin[Qj(/ — s)]dl (19) 

rt 


(nS)i' 


- s)fi(s)ds 


( 20 ) 


with K 0 (t — l) = sin[fi(t — /)], f 0 (t) = fo(t), and fi(s) given by Eq. (ITKlh Hence, 
Vn < N, Eq. (1T5li for x(t) can be rewritten as follows: 


x{t) = f n (t ) + ( n 


i=l \i =o 


A \ A-i 

A j ~dT 


Ki(t - s)X i _ 1 (s)ds 


+ ^ J K n (t - s)X n+l (s)ds . 


( 21 ) 
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Proof: The proof is by induction. First of all we show that for n — 1, Eq. m is 
correct. Substituting Eq. (fTBli for Xi(t) in Eq. (fTT!) one finds 

x(t) = f 0 {t ) + sin[fi(t - s)]/i(s)ds 

D 2 A DD f t 

+ ^ s ) x ( s ) ds + qq 1 / s)A 2 (s)ds. (22) 

i J o i *y o 

This equation can be easily recast in a form like Eq. (OTT) . Assume now that Eq. (EH 
is true for a generic n < N, and substitute Eq. ([Til for X n+1 (s) in the second line 
of Eq. EB . After some simple manipulation one can show that the Eq. (T2TT) reads 
now as follows: 


«M = jUM + g(ri 

/ydi n \ A 

+ I ]f[ ^ ) / ^ n +\{t ~ s )X n +2 (s)ds 
\i =o 4 v “'° 


s)AVi(s)ds 


(23) 


that is exactly Eq. (T2Tk for n —> n + 1. This completes the proof. □ 


It is important to note that equation EH is exact. Although only the equations 
for the first n have been substituted into x(t ), the term depending on X n+ \ of 
Eq. EB encodes the information regarding the evolution of the remaining N — n 
modes. It is easy to prove that, if one substitutes the equations for all the N AT, 
the following Corollary holds true: 


Corollary: Let the hypothesis of Theorem 1 be true, and let n = N. Hence, 


N 


X{t) = fN(t)+^2 If 


i= 1 V=0 


A \ A-i 

A j “dT 


I<i(t - s)Xi_i(s)ds 


(24) 


Proof: By definition, D^ = 0. Accordingly, the second line of Eq. EB is null. □ 

This equation determines how the dynamics of x(t) is affected by the full set of 
Xf f]y(t) accounts for the initial conditions of the collective modes, while the second 
term is a purely non-Markovian contribution which involves the whole past evolution 
of the collective modes. It is important to underline that each AT contributes to x(t) 
via K i+ i. This is a crucial feature of the dynamics that will play a fundamental role 
in the following calculations. 
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In Eqs. (I2T|) and (I24p the dependence on x(t) is not completely explicit. Recalling 
that X 0 = x we rewrite such equations as follows: 

x (t) = Typr [ Ki(t - s)x(s)ds + F N (t) (25) 

\i\li J 0 

where the function F N (t) collecting all terms that do not depend on x is defined as 
follows: 

Ki(t - s)Xi-i(s)ds 
( s)ds (26) 

Ki(t - s)X i ^ 1 (s)ds , (27) 

where the definitions come from Eq. (I2T1) and Eq. (OH) respectively. Note that if 
one chooses n — N in the first equation, one recovers the second one as expected. 
Eq. 051) explicitly shows that the dynamics of x is ruled by an integral equation. 
Since the kernel K\(t — s) is a linear combination of two sine functions, Eq. 051) can 
be solved using standard techniques^. The solution reads 

D 2 f* 

x(t) = F N (t ) H- 7-2 - 27 / (A* 2 sin[^i (t - s)] - /ii sin [jj. 2 (t - s)]) F N (s)ds 

hlh2l/ i 2 — hlJ JO 

(28) 

where 

Mi ,2 = J\ (^ 2 + ± \/A) , A = (Q 2 - Q 2 ) 2 - 4D 2 (29) 

In order to avoid multi-valued x(t), n 1,2 have to be real, i.e the condition A > 0 has 
to hold true. Equation OHl) is the exact solution of our problem, as it displays the 
dynamics of an harmonic oscillator under the influence of a chain of N harmonic 
oscillators. One can identify two different type of contributions: a diffusive one, 
given by those terms contained in F N (t) that depend on the initial conditions qj.{ 0); 
and a purely non-Markovian one, given by the integral terms, and which depends on 
the interaction among the chain modes. This equation gives a description of non- 
Markovian dynamics which is equivalent to the Generalized Langevin Equation^. 
Furthermore, Eq. (1281) plays an important role because it has the added value of 
allowing for an effective description of non-Markovian dynamics, i.e. a description 
in terms of a smaller number of degrees of freedom. 


N 


ivw=/„<*)+£in#' A_i 


i =2 \i =0 


0 / D 


I JO 


n £)/* **<*-)*- 


N 




i= 2 \l =0 


a / d 


I JO 
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III. EFFECTIVE NON-MARKOVIAN DYNAMICS 


In the study of open quantum system the number N of constituents of the envi¬ 
ronment is typically assumed to be very large. Controlling the dynamics of so many 
constituents is not possible. One then has to consider approximated dynamics, or 
in many practical applications one needs to exploit numerical methods. Aim of this 
section is to understand whether it is possible to describe with a good approxima¬ 
tion the full dynamics of an open system, by a particle interacting with a smaller 
bath (in our case, a shorter chain). In fact, we will prove that there always exist a 
time-scale such that the true dynamics x(t ) of the system is well approximated by 
truncating the original chain after n < N oscillators. Let x (n) (t) denote the the ex¬ 
act evolution of x interacting with a chain of n oscillators. Observe that truncating 
the chain after the n-th oscillator corresponds to setting D n = 0. Accordingly, the 
evolution of x (n) (t) is given by Eq. (OH1) with N replaced by n. It proves useful for 
the forthcoming discussion to define the following function: 


ei(n,t) ■= F N (t ) - F n (t) 


(30) 


Comparing Eq. (|26li for ibv(f), and Eq. (|27li for F n (t), one can easily show that 


•■■-'-(StK 


X n (t 


(31) 


We define a further function that measures the error made by considering the 
truncated dynamics instead of the full one: 


e(n,t) := \x(t) - x (n) (t)\ (32) 

Exploiting Eq. (OSli one Ends that 

e(n , t ) = | ei (n, t ) + e 2 ( n , t) \ (33) 


where 

e 2 (n,t) = 


D 2 


(/i 2 sin[/ii(f - s)] - sin[/i 2 (t - s)]) ei(n, s)ds (34) 


hiMh! ~ vl) Jo 

Our aim is to obtain an upper bound on the error function e(?r, t). Since e(n,t) 
strongly depends on the features of the kernels K n (t — s ), in order to reach our goal 
we need to understand their features. 
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A. Kernels structure 


From the definition of Eq. ([TUI) , one easily understands that each AT consists of i 
nested integrals of sine functions. This structure, that is due to the harmonic feature 
of the chain, turns out to be crucial for the analysis of the system. Indeed, the h-th 
derivatives of AT have the following remarkable properties. 


Theorem 2: Let K,(t — s) : M 2 —» M be the kernels defined in Eq. (THJ1) . and let 
Kf ] (t — s) be their k-\\i derivatives with respect to t. Hence, the following equations 
hold true: 


K\ k \t — s) = J K^[(t — l) sin[f 2 j(/ — s)]dl Vfc < 2i — 1 

k 


-(k) 


(35) 


- s) = sm[Sij(i - s)] ^(-l) ,+ *JCi?i 1) (0)U 

3 = 1 
rt 

q2i+i) 


2k—2j 


+ / K>_i l, {t-l)sm[Qi{l-s)]dl Vk>i 


(36) 


fc -1 


K 


(2fe—1) 


(t- s )= CO S [ft.(« - S )] ^(-l)> + '=- 1 Afy 1| (0)fif- 1 -^ 


3 =1 


+ j KfJ+V(t - l) sm[n t (l - s)]dl Vk>i + 1 
Proof: The proof is by induction. We recall the definition of AT(t — s): 


(37) 


Ki(t — s) = / A' 0 (t — /) sin[fli(/ — s)]dl = / sin[H(t — /)] sin[fli(Z — s)]dZ (38) 

J S J S 

Differentiating this equation it is easy to check that k[ 1 ^, k[ 2> and satisfy 
the system above. The proof for higher derivatives can easily be done differentiating 
recursively. Assume now that the system (j35[) - (j3TI) holds true, and note that this im¬ 
plies that Kf :) (0) = 0 for all k < 2i — 1 . Accordingly, if one iteratively differentiates 
the definition (fT9l) for AT+i, one finds 


Ki+i(t — s) = I K\ K, (t — l ) sin[Dj(/ — s)]dl Wk < 2i + 1 


"(k) 


(39) 


that is Eq. (1351) for K i+i . Differentiating this equation one can easily prove that 
Eqs. (j36|) - ([3T|) hold true for K l+ \ as well. □ 

This Theorem shows some interesting features of the kernels derivatives which 
will prove essential. First of all, Kf ] (0) = 0 not only for all k < 2i — 1, but also for 
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all even k. Moreover, iterating Eq. (HT71) one can show that 


0 ) 



E 

aj=0 

T, OLj =k—l—i 



(40) 


We exploit Theorem 2 to nnderstand the time behavior of the kernels, by ex¬ 
panding I\i in Taylor series. First of all note that, since the first 2 i and all the even 
derivatives of K{ are null, the kernels Taylor series start with the index i — 1, and 
they display only odd terms. 


K i(t - s) = Yl 

k =0 


Af’( 0) 

k\ 


{t - S ) fc 


OO 


E 

k=i +1 


K\ 2k ~ l \ 0) 
{2k — 1)! 


{t-s) 


2k-l 


(41) 


This equation gives an interesting insight on the evolution of x{t ) and on how the 
non-Markovian behavior emerges. We recall Eq. (OTli which shows that each A* 
influences x{t) via K i+i . Accordingly, if we consider t < 1, the smaller i the earlier 
X t contributes to the dynamics of x. At the beginning, only Xi gives a relevant 
contribution to the non-Markovian term of the dynamics of x{t ), i.e. the integral 
term of Eq. (|24l) . and as time grows, also other modes enter into the game. In other 
words, for short time scales only the first oscillators of the chain contribute to the 
non-Markovian dynamics. 


B. Short-time approximation 

From the discussion above one can also infer another fundamental property of 
the system: if the time scale is short enough, further oscillators can be neglected 
since their contribution is very small. Accordingly, one can effectively describe the 
dynamics of the full system by a truncated chain. The time scale, the number 
of oscillators in the truncated chain, and the error one makes using the effective 
dynamics, are strictly connected quantities. The next Theorem proves how these 
quantities are related. 

Theorem 3: The error function e{n,t) defined in Eqs. (I30l) - (l34l) is bounded from 
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above as follows: 


N 


k =1 


(n,t) < ^ | PnM\ ( |?*(0)| + ) t 2n+2 


|4(0)| \ 

Wfc / 


cosh (VE’L o) D2f4 cosh (+ EILo 


(2n + 2)! 


(2n + 6)! 


(42) 


Proof: We apply the triangular inequality to Eq. (H31h and we evaluate indepen¬ 
dently the two contributions of |ei(n, t)\ and \e 2 (n,t)\. The first one reads 


..■"dfi 8)1 

f n \ 1 

nw 1 ei^)i(w°)k 


K n {t ~ S)| • |X n+ i(s)| ds 

l4(o)| 


yl = 0 


k =1 




| K n (t-s)\ds (43) 


where the second line is obtained by expressing X n+1 in terms of the q k by means of 
Eqs. © and (jSlh Since the q k are independent, they evolve with linear combinations 
of sines and cosines which have been bounded by 1. One can then exploit Eq. 
and obtain 

rt °° Kn jl \o) 

2jT 

j=n +1 


rt ^ 

/ \K n (t-s)\ds< V] 
Jo 




OlPd E t2 ‘ 


n \ j - n ~ 1 

l2 


< 


J =0 


n* 


j=n +1 

f.2n+2 


2 j! 


Eh 


. i=1 


cosh t. 


Eh 


(44) 


, (2« + 2)! V N"=t 

The second inequality is obtained by substituting Eq. (HU for Kn J 1 ' ) (0) and by 
observing that this equation is essentially (Xq=owith all the coefficients 
set to 1. The third line comes from a change of variable on the sum index and 
extending the series to zero. Substituting this result in Eq. (1421) one finds 

mw)jzz cosh u 


hi 


< ^2\P n (u k )\ I \q k {0)\ + 


k =1 


uj k J (2 n + 2)! 


\ 


Eh 


(45) 


i=l 


The contribution of \e 2 (n,t)\ to Eq. (|42li is obtained as follows. First of all one 
observes that the integral kernel of e 2 (n,t) in (T34l) is of the same order as K\\ 
/i 2 sin[/ii(f - s)] - /ii sin- s)] 


(4 ~ at 


:= K\(t-s). 


(46) 
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Substituting the definition of e±(n,t) in Eq. (1771) . one finds that X n+ \ contributes 
via a kernel K n+2 of the order n + 2. The final result is obtained following the same 
procedure as for |ei(n, t)\. □ 

As expected Eq. (1721) displays a dependence on the initial conditions of the bath 
position operators. Since these quantities are in principle not known, to overcome 
this issue one traces over the bath degrees of freedom. We assume that the bath 
initial state is in thermal equilibrium at temperature T : 



(47) 


where Z guarantees that Tr p = 1, and (3 = (ksT) 1 with ks the Boltzmann constant. 
We prove the following Theorem: 


Theorem 4: If the initial state of the bath is given by Eq. (1471) . the error is in 
average bounded from above as follows: 


<e(M)> < 



^2n+2 _ 


cosh L> 2 £ 4 cosh [t^Vt 2 + Of + ]T" =0 ttfj 

(2n + 2)! + (2n + 6)! 


(48) 


Proof: We consider the term in Eq. (1T71) that depends on the initial operators, and 
we average it over the initial state (1471) : 



One can easily perform the two X-dimensional Gaussian integrals and obtain 


N 


8 k B T \P n (uk)\ 


IT 


E 

k= 1 




(50) 


Replacing this expression in Eq. (1771) one finally obtains Eq. (I48[) . □ 

Theorem 4 shows how e, t and n relate to each other at a given temperature. 
This is a fundamental result as it explains in which sense and under which limits 
the truncated dynamics can be considered an effective dynamics for the full system. 
Moreover, the result is versatile and can be easily adapted to different uses, since one 
can fix two of the mentioned quantities in order to obtain the third one. For example, 
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one might be interested to know how big is the error after a time t, truncating the 
chain to n oscillators. On the other side, one might choose to fix a value for the 
error at a given time, and derive how many oscillators are needed in order to satisfy 
such condition. It is not easy to invert analytically Eq. fl48l) . but it is rather easy to 
evaluate it numerically. 

IV. CONCLUSIONS 

In order to derive an effective description of non-Markovian open systems dy¬ 
namics, we considered a chain representation of the environment. After deriving 
analytically the exact dynamics of the problem, we proved the existence of an effec¬ 
tive description for such a dynamics and we provided its explicit form. The proof is 
based on a peculiar feature of the integral kernels of the dynamics, deeply related to 
the harmonic feature of the model. In particular, this intrinsic structure shows that 
“far” oscillators can be neglected if the timescale of the dynamics is short enough. 
Furthermore, we gave an upper bound of the error one makes when considering 
such an approximated dynamics instead of the exact one. Such bound depends on 
the number of oscillators in the truncated chain, on the timescale and on the tem¬ 
perature of the system, setting a strong relationship among these quantities. This 
achievement gives a strong basis for future investigations on non-Markovian systems 
both at the analytical and numerical levels. For example, our result provides a way 
to determine the number of modes in a chain necessary to keep the error small, or 
the timescale at which a numerical simulation can be considered accurate. 
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